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^■f^ ' We study algorithmic aspects of bending wires and sheet metal into 

a specified structure. Problems of this type are closely related to the 
question of deciding whether a simple non-self-intersecting wire structure 
(a carpenter's ruler) can be straightened, a problem that was open for 
, several years and has only recently been solved in the affirmative. 

• ' If we impose some of the constraints that are imposed by the manufac- 

' turing process, we obtain quite different results. In particular, we study 

the variant of the carpenter's ruler problem in which there is a restriction 
that only one joint can be modified at a time. For a linkage that does 
not self-intersect or self-touch, the recent results of Connelly et al. and 
Streinu imply that it can always be straightened, modifying one joint at 
a time. However, we show that for a linkage with even a single vertex 
^■—^ , degeneracy, it becomes NP-hard to decide if it can be straightened while 

■ altering only one joint at a time. If we add the restriction that each joint 

' can be altered at most once, we show that the problem is NP-complete 

^SJ . even without vertex degeneracies. 

• In the special case, arising in wire forming manufacturing, that each 

"'^ I joint can be altered at most once, and must be done sequentially from one 

^ . or both ends of the linkage, we give an efficient algorithm to determine if 

a linkage can be straightened. 
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1 Introduction 



The following is an algorithmic problem that arises in the study of the man- 
ufacturability of sheet metal parts: Given a flat piece, F, of sheet metal (or 
cardboard, or other hendable stiff sheet material), can a desired final polyhe- 
dral part, P , be made from it? The 2-dimensional version is the wire-bending 
("paperclip") problem: Given a straight piece, F, of wire, can a desired simple 
polygonal chain, P, be made from it? This problem also arises in the fabrica- 
tion of hydraulic tubes, e.g., in airplane manufacturing.^ In both versions of the 
problem, we require that any intermediate configuration during the manufac- 
ture of the part be feasible, meaning that it should not be self-intersecting. In 
particular, the paperclips that we manufacture are not allowed to be "pretzels" 
- we assume that the wire must stay within the plane, and not cross over itself. 
See Figure |l| for an illustration. We acknowledge that some real paperclips are 
designed to cross over themselves, such as the butterfly style of clip shown in 
the figure. 

Our problem is one of automated process planning: Determine a sequence 
(if one exists) for performing the bend operations in sheet metal manufacturing. 
We take a somewhat idealized approach in this paper, in that we do not attempt 
to model here the important aspects of tool setup, grasp positions, robot motion 
plans, or specific sheet metal material properties which may affect the process. 
Instead, we focus on the precise algorithmic problem of determining a sequence 
for bend operations, on a given sheet of material with given bend lines, assuming 
that the only constraint to performing a bend along a given bend line is whether 
or not the structure intersects itself at any time during the bend operation. 

Note that the problem of determining if a bend sequence exists that allows a 
structure to unfold is equivalent to that of determining if a bend sequence exists 
that allows one to fold a fiat (or straight) input into the desired final structure: 
the bending operations can simply be reversed. For the remainder of the paper, 
we will speak only of unfolding or straightening. 

1.1 Motivation and Related Work 

Our foldability problem is motivated from process planning in manufacturing 
of structures from wire, tubing, sheet metal, and cardboard. The CAD/CAM 
scientific community has studied extensively the problem of manufacturability 
of sheet metal structures; see the thesis of Wang |34| for a survey. Systems have 
been built (e.g., PART-S ||l2| and BendCad [|7|) to do computer-aided process 
planning in the context of sheet metal manufacturing; see also ||, HI ^ . 

See for a motion planning approach to the problem of computing folding 
sequences for folding three-dimensional cardboard cartons. Considerable effort 
has gone into the design of good heuristics for determining a bend sequence; 
however, the known algorithms are based on heuristic search (e.g., A*) in large 

^We thank Karcl Zikan for introducing to us the hydrauhc tube bending problem at Boe- 
ing's factory. 
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state spaces; they are known to be worst-case exponential. (Wang Q cites the 
known complexity as 0(n!2").) 

Our work is also motivated by the mathematical study of origami, which 
has received considerable attention in recent years. In mathematics of origami, 
Bern and Hayes ||] have studied the algorithmic complexity of deciding if a given 
crease pattern can be folded flat; they give an NP-hardness proof. Lang 
^ gives algorithms for computing crease patterns in order to achieve desired 
shapes in three dimensions. Other work on computational origami includes 
|l|, 1^, H, 111. A closely related problem is that of flat foldings of 

polyhedra. It is a classic open question whether or not every convex polytope in 
three dimensions can be cut open along its edges so that it unfolds flat, without 
overlaps. Other variants and special cases have been studied; see [|[ ^, ||, p5|, 

Finally, we are motivated by the study of linkage problems; in fact, in the 
time since this paper was first drafted, the carpenter's ruler conjecture has 
been resolved by Connelly, Demaine, and Rote and Streinu [|^: Any 
(strongly) simple polygonal linkage with fixed length links and hinged joints, 
can be straightened while maintaining strong simplicity (i.e., without the link- 
age crossing or touching itself). (They also show related facts about linkage 
systems, e.g., that any simple polygonal linkage can be convexified.) In fact, 
Streinu | |3^ gives an algorithmic solution that bounds the complexity of the 
unfolding and is somewhat more general than the slightly earlier results of . 
These results imply that any (strongly simple) paperclip can be manufactured if 
one has a machine that can perform a sufficiently rich set of bending operations. 
For a recent overview of folding and unfolding, see the thesis of Demaine |p3|. 
Earher and related work on hnkages includes |0, ^, |2^, 23, ^ ^ Oi 



hardness results are particularly interesting and relevant in light of these new 
developments, since we show that even slight changes in the assumptions about 
the model or the allowed input results in linkages that cannot be straightened, 
and it is NP-hard to decide if they can be straightened. 



1.2 Summary of Results 

(1) We show that it is (weakly) NP-complete to determine if a given rectilinear 

polygonal linkage can be straightened, under the restriction that only one 
joint at a time is altered and each joint can be altered only once (so the 
joint must be straightened in a single bend operation). A consequence is 
that the more general sheet metal bending problem is hard as well, even 
in the case of parallel bend lines and an orthohedral structure P. 

(2) We prove that it is (weakly) NP-hard to determine if a given polygonal 

linkage can be straightened if there is a vertex degeneracy, in which two 
vertices coincide. Here we again assume that only one joint can be altered 
at a time, but we do not assume that a joint is altered only once, so we 
may make any number of bends at any particular joint. 

(3) We give efficient algorithms for determining if a given bend sequence is 

feasible, assuming only one joint is altered at a time, and for determining 
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if certain special classes of bend sequences are feasible. In particular, we 
give an efficient (O(nlog^n)) algorithm for determining if a polygonal 
linkage can be straightened using a sequential strategy, in which the joints 
are completely straightened, one by one, in order along the linkage. We 
also give efficient polynomial-time algorithms for deciding whether there 
is a feasible bend sequence that straightens joints in an order "inwards" 
from both ends or "outwards" towards both ends. (Such constrained bend 
sequences may be required for automated wire-bending machines.) These 
results will be made more precise in Section |[ 



2 Preliminaries 

The input to our problem is a simple polygonal chain (linkage), P, with vertex 
sequence [bo, 61, 627 • ■ • > ^n+i)- The points bo and 6„+i are the endpoints of the 
chain, and the n vertices 61, . . . , 6^ are the bends (or joints). The line segments 
bibi+i are the edges (or links) of P. The edge bibi+i is a closed line segment; 
i.e., it includes its endpoints. We consider the chain P to be oriented from 60 
to bn+i, and we consider each edge of P to have a left and a right side. Each 
bend bi has an associated bend angle 9i S (0, 27r], measured between the right 
sides of the two edges incident on bi. 

The chain P is strongly simple if any two edges, bibi+i and bjbj+i, of P that 
are not adjacent {i ^ j) are disjoint and any two adjacent edges share only 
their one common endpoint. We say that P is simple if it is not self-crossing 
but it possibly is self-touching, with a joint falling exactly on a non-incident 
edge or another joint; i.e., P is simple if it is strongly simple or an infinitesimal 
perturbation of it is strongly simple. 

We consider the chain P to be a structure consisting of rigid rods as edges, 
whose lengths cannot change, connected by hinged joints. When a bend oper- 
ation is performed at joint bi, the bend angle 9i is changed. Throughout this 
paper, we assume that the only bend operations allowed are single-joint bends, 
in which only one bend angle is altered at a time. We establish the convention 
that when a bend operation occurs at 6^, the subchain containing the endpoint 
bo remains fixed in the plane, while the subchain containing bn+i rotates about 
the joint bi. This convention allows us to have a unique embedding of a partially 
or fully straightened chain in the plane. 

A bend operation is complete if, at the end of the operation, the bend angle 
is tt; we then say that the joint has been straightened. A bend operation that 
is not complete is called a partial bend. A sequence of bend operations is said 
to be monotonic if no bend operation increases the absolute deviation from 
straightness, \9i — tt\, for a joint bi. If all joints of P have been straightened, 
the resulting chain is a straight line segment, F, of length J2^=o \bibi+i\, where 
\bibi+i \ denotes the Euclidean length of segment bibi+i. By our bend operation 
convention, one endpoint of F is bo, and F contains the segment bobi (which 
never moves during bend operations). 

For S* C B, we let P{S), denote the partially straightened polygonal chain 
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having each of the bends bi £ S straightened (to bend angle tt), while each of 
the other bends, bi (ji S, is at its original bend angle 9i. Thus, in this notation 
P{B) = F and P(0) = P. We let P{S;i,0), for 1 < i < n with i ^ S, denote 
the chain in which each bend bj £ 5 is at bend angle tt, bend bi is at angle 0, 
and all other bends bj ^ 5* are at their original bend angles dj. We say that 
chain P{S) or P{S; i, 0) is feasible if it is a simple chain. 

We say that bend bi is foldable (or is a feasible fold) for P{S) if P{S; i, 9) is 
feasible for all 6 in the range between tt and 9i (more precisely, for all 9i < < tt, 
if 9i < TT, or for all tt < 9 < 9i, if Oi > tt). If bi is foldable, then it is possible 
to make a complete bend at bi, meaning that the joint can be straightened in 
a single operation without causing the chain to self-intersect. We say that a 
permutation a = (ii, 22, . . . , in) of the indices {1, 2, . . . , 71} is foldable for P if, 
for j — 1, 2, . . . , n, joint bi- is foldable for P{{bi^, . . . , bi._-^}), i.e., if P can be 
unfolded into the straight segment F using the bend sequence a (so that, by 
reversing the operations, P can be manufactured from F using the reverse of 
the bend sequence). 

The Wire Bend Sequencing problem can be formally stated as: Deter- 
mine a foldable permutation a, if one exists, for a given chain P. 

This paper studies the Wire Bend Sequencing problem for polygonal 
chains in the plane. We note, however, that our results have some immediate 
imphcations for the Sheet Metal Bend Sequencing problem, which is de- 
fined analogously for a polyhedral surface P having a pattern B of bend lines 
(creases), each of which must be straightened in order to flatten P into a flat 
polygon F. Specifically, the hardness of the Sheet Metal Bend Sequencing 
follows from the hardness of the Wire Bend Sequencing, which can be seen 
as a special case of the sheet metal problem in which F is a rectangle and the 
bend lines B are all segments parallel to two of the sides of F and extending all 
the way across F. 

We give an example in Figure |^ of some common paperclip shapes, (a)-(c). 
We also show an example, (d), of a 5-link paperclip that cannot be straight- 
ened using complete bends, for any permutation a of the bends. Finally, we 
show an example of a 6-link paperclip for which the foldable permutations are 
{(1, 5, 4, 3, 2), (1, 5, 4, 2, 3)}; we show the sequence of bends, with the intermedi- 
ate structures, for the permutation a = (1, 5, 4, 3, 2). 
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Figure 1: Examples of paperclips: (a) and (b) are standard versions, which 
are readily straightened, (c) is a "butterfly" paperclip, which is not a planar 
structure and is not among the wire structures considered in our two-dimensional 
model, (d) shows a 5-link paperclip that cannot be straightened using complete 
bends in the plane, (e) shows a 6-link structure that can be straightened, 
e.g., using the bend sequence animated below it for the bend sequence a = 
(1,5,4,3,2). 
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3 Hardness Results 



3.1 Complete Bends 

Our first result shows that if we require bends to be complete, as in our spec- 
ification of the Wire Bend Sequencing problem, the problem of deciding if 
there is a feasible bend sequence is NP-complete. 

Theorem 3.1 Wire Bend Sequencing is (weakly) NP-complete, even if P 
is rectilinear. 

Proof: We prove NP-completeness, even in the case that we are restricted 
to a special class of bend sequences, namely, those that can be written as the 
concatenation of up to four monotone subsequences of the index set {1, . . . , n}. 
Below, we refer to each subsequence of bends as a monotone pass over the chain, 
going from one end to the other, performing a specified subset of complete bends. 

Our reduction is from Partition: Given a set S* of n integers, a^, which sum 
to A — ^ j Qi , determine if there exists a partition of the set into two subsets 
each of which sums to A/ 2. 



partition 
elements 




Figure 2: Proving hardness of the Wire Bend Sequencing problem for recti- 
linear chains: Frame and key. 



The key idea of our construction uses two components, as shown in Figure ||: 
One is a rigid "frame" that can only be unfolded if one end of the chain can be 
removed from within this frame. The other component is a "key" that encodes 
the partition instance. Collapsing the key is possible if and only if there is 
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a partition of the integers into two sets of equal sum. The total number of 
segments will he £ — 26 + 4n; we write bi (i = 0, 1, . . . , 26 + 4n) for the vertices, 
and Si — for the segments. For any point in time, we refer to the 

position of a joint bi by its coordinates {xi,yi). When discussing some of the 
relative distances, we use dooibi, bj) = ma,x{\xi — Xj\, \yi — yj\}- 

More precisely, the frame consist of 13 segments, si — {bo,bi), . . . , S13 = 
{bi2, bis), as shown in the figure. Segment lengths are chosen such that the size of 
the frame is Q{L), with minimal coordinate differences doc{ba, ^13), doo{bi, ^'12), 
rfoo(&2,&ii), docibzMo), doa{bi,bg), doo{bhM), docib&M) being e(e), where e = 
l/{n^L^). The "key" consists of 13 + 4n segments, S14 — (613, fei4), . . . , S26+4n = 
(J'25+4n,^26+4n)- For I = 0, . . . , 4n + 2, the "auxiliary" segments Si5+4i, Si6+4i, 
Si7+4i have length e, while the "partition" segments Si8+4i have length a^. The 
long "positioning" segments S14, S15, S25+4m S26+4n have lengths L, L/3, L/3, 
and L — 'ine, respectively; they guarantee that the partition segments must have 
a particular relative position when removing the key. We choose the scale to 
be such that L/A > A, for technical reasons that will become clear later in 
the proof. As indicated in the figure, the initial position of each key segment 
Si, i = 14, . . . , 26 + 4n has a;-coordinate 0:13 or X13 +e, with a horizontal distance 
of = 2e from S4. Moreover, 614 is positioned at a vertical 

distance of 2/14 — — ri^e = l/{nL^) above 64. 

The purpose of the auxiliary segments is as follows. As shown in Figure ^ 
we have two types of joints in the figure: the "ordinary" ones (indicated by solid 
black dots) form the frame and can only be accessed once. The "quadruple" 
ones (indicated by hollow dots in Figure ||) consist of the four simple joints at 
three consecutive auxiliary segments; they are found along the key as described. 
These quadruple joints make it possible to simulate opening and closing such a 
joint a limited number of times. 

Now assume that there is a partition S — Si L) S2, such that X^ieSi ~ 
SieS2 order to see that the key can be removed from the frame we first 

convert it into the "stair" configuration shown in Figure |^: We make one mono- 
tone pass over the chain towards the key end, and straighten one ordinary joint 
per quadruple joint whenever this joint separates two segments from different 
Si. Thus, segments corresponding to numbers in Si will be horizontal, while 
those for numbers in ^2 will be vertical. In order to keep the number of mono- 
tone passes limited to four, during this first pass we also straighten two ordinary 
joints per quadruple joint separating two segments from the same Si, as shown 
in Figure ^(b). 

Making a similar monotone pass, we can convert the stair into a "flat har- 
monica", as shown in Figure ^, with segments from 5*2 pointing "down", i.e., 
yi < yi-i, and segments from pointing "up", i.e., yi > yi^i- By as- 
sumption about the partition, the positions of endpoints 613 and 5i8+4n satisfy 
rfoo(fei8 7 fei8+4n) < ine, and 2L/3 — 3e < 1/13 — yis < 2L/3-|-3e, i.e., both 618 and 
6i8+4n are roughly 2L/3 below 613. Altogether, the position of the last segment 
S26+4n of length L in the chain will differ by at most 0{n£) from the verti- 
cal position of segment S14, with all other segments strictly in-between. This 
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(a) 



(b) 



Figure 3: Turning tlic key into a stair: (a) An intermediate stage of the mono- 
tone pass, (b) The stair configuration at the end of the monotone pass, with 
details of the state of quadruple joints. 




n'£±O(n£) = 0(l/nL)«l 



Figure 4: Turning the stair into a harmonica of small width and length L. (The 
horizontal width is not drawn to scale in order to show details.) 



9 



collapsed structure can be rotated about 613 without colliding with any frame 
segments. Then it is easy to open up the remaining frame (by straightening 612, 
611, 610, bs, &7, be, 65, 64, 63, 62, bi as one monotone pass, skipping 6g.) Finally, 
the resulting monotone chain can be straightened in one last monotone pass. 

Conversely, assume now that the chain can be straightened. See Figure ^j. 
It is clear that 613 must be straightened before any other joint in the set 
{61, . . . , 612}. In order to avoid hitting vertex 64 during this motion, any part of 
the key to the right and below &13 must be strictly within the circle C of radius 
r = + l/(nL2))2 ^ ^2 < 7^ + 2/{nL'^) around 613, where r is the distance 
between &13 and 64 (see Figure |[) The following technical arguments show 
that at this time, segment S26+4ri has to be in a vertical position that basically 
coincides with su, which is only possible in case of a feasible partition. 




critical segment s 26+4n 



'12 
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'11 



Figure 5: When straightening joint 613, the key must be fully contained in the 
shaded circle of radius r < L + 2/[nL?). This forces a particular position of 
segment S26+4n- 



When starting the rotation about 613, S26+4n is an axis-parallel segment of 
length L — Sne > L — l/{nL'^). The rigid frame and the closeness of 614 and 64 
ensure that segment S26+4n cannot lie to the left of S14, implying that S26+4n 
can only lie within the quarter circle of radius r below and to the right of 613 
when 613 is straightened. 

Let 6min be one of the two points in {625+4„, &26+4n} that is not further from 
613 than the other, and let ^max be the other point. If the vertical distance j/13 — 
2/min is greater than i/7/ {nL), it follows that the Euclidean distance between 



10 




and &13 is at least 

a contradiction to the assumption that S26+4n is fuUy contained in C . Now, 
using our assumption that I//4 > ^, we know that 625+4T1 and &14 are connected 
by a polygonal chain of length strictly less than L/3 + L/4 + L/3 = llL/12, 
implying that 625+4n has Euclidean distance at least L/12 from 613, so 626+4ri = 
6niin and 625+4ri = ^^max- As fe26+4n IS within yJYJ{nIj) of 613, it follows that 
^26+4n has Euclidean distance at least L — a/7/ inV) from 614. If S26+4ri were 
horizontal, then the Euclidean distance between &25+4n and 614 would be at 

least ^ (^L - / {nL)^ + (i - l/(ni2))2 > llL/12, a contradiction. Hence, 

■S26+4n must bc vertical. Just as we derived for the vertical distance between 613 
and 626+4n, it follows for the horizontal distance that a;26+4n~2;i3 < \/7/ {nL) <C 
1. 

Now observe that when starting the rotation about 613, all partition segments 
must be strictly between su and the narrow strip between su and S26+4ti, 
meaning that they are all vertical. Let Si be the set of "upwards" partition 
segments Si with > Ui-i, and 5*2 be the set of "downwards" partition segments 
Si with Hi < Ui^i. As |y24+4n - 2/15 1 = Q{ne) and |j/25+4n - Vul = Q{n'^e), we 
conclude that the integral total length of upwards segments equals the integral 
total length of downwards segments. 

This means that X^ieSi ^« ~ X]iGS2 have a feasible partition. This 

completes the proof. □ 



3.2 Partial Bends 

Now we consider the case in which each joint may be changed an arbitrary 
number of times during the straightening operations, while still making single- 
joint bends (bending only one joint at a time). This version of the problem 
is closely related to the carpenter's ruler problem studied by |l^, |3^. In the 
context of our study on folding, there may be the additional requirement of 
using only monotonic bend operations, e.g., to avoid work- hardening the wire, 
possibly causing it to break. We begin with the following observation about 
the sufficiency of monotonic single-joint bends; see also the discussion on p. 9 
of Demaine's thesis [|l3| . 

Theorem 3.2 Any strongly simple polygonal chain P can he straightened using 
a finite number of monotonic single-joint bends. 

Proof: Consider the set S of points in n-dimensional joint-angle space that 
correspond to strongly simple embeddings of the linkage. A single-joint bend 
corresponds to axis-parallel motion in joint-angle space. If self-touching is pro- 
hibited, S is an open set; note too that S is bounded. By Streinu's result pa. 
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there is an opening motion of the chain that consists of a finite number of in- 
dividual monotonic moves. Such an opening motion corresponds to a path, 11, 
in S, comprised of a finite number of arcs, each corresponding to a monotonic 
move. Let e be the Euchdean distance between path 11 and the boundary of S; 
since S is open, we know that e > 0. Then we can replace each arc of the path 11 
with a finite sequence of axis-parallel moves of size e/2, yielding a straightening 
that uses single-joint bends. □ 

We will refer to a sequence of small individual moves that mimics an overall 
large-scale motion of several joints as "wiggly" , since the overall motion may be 
achieved through back-and-forth motions of individual segments that gradually 
change individual angles. 

The following results show that allowing even a single point of self-incidence 
along the linkage changes the overall situation quite drastically. 

Lemma 3.3 There are polygonal chains P with a single vertex-to-vertex inci- 
dence that cannot be straightened using partial single-joint bends. 

Proof: See Figure ^ The chain has eight joints (labeled 6o, . . . ,67) and 
seven segments (of the form Si — (bi-i,bi)). The endpoint bo coincides with 
joint 65. It is easily checked that none of the joints 5i, . . . ,64 can be changed 
without causing a self-intersection: Assume that there is a feasible motion of 
a joint bi with < i < 5. Then the points bo and 65 would move away from 
each other along a circle around bi. Without loss of generality, assume that 
65 remains in place, while bo is moving. Now consider the first such rotation 
that starts with 60 and 65 coinciding, and that avoids a crossing of si with both 
S5 and sq. If bo moves clockwise around 6^, it is easy to see that the angle 
between {bo,bi) and S5 must be at least it/2 when starting the motion, or else 
si and S5 intersect. If 60 moves counterclockwise around bi, the same follows 
for the angle between {bo,bi) and sg. Therefore, the center of rotation must 
lie within the shaded region shown in the figure. (The cone to the left of 65 is 
feasible for clockwise rotation, while the cone to the right of 65 is feasible for 
counterclockwise roation.) However, none of the joints 5i, . . . , 64 lies inside of 
this feasible region. It follows that bo, ... ^b^ form a rigid frame, as long as the 
angle at 65 stays smaller than 7r/2. 

On the other hand, it is easy to see that bj cannot be removed from the 
pocket formed by 61, 62, and 63 if only the two remaining "free" joints 65 and 
be can be changed. The claim follows. □ 

If 60 and 65 have some positive distance, then the frame can be opened along 
the lines of the approach in or |Q by gradually straightening 65, be, bi, 62, 
and 63, so that the "zig-zagging" part between bo and 63 pushes left, while be 
swings around 65. 

Using the frame as a gadget, we can show the following: 

Theorem 3.4 It is NP-hard to decide if a polygonal chain P with a single 
vertex-to-vertex incidence can be straightened by arbitrary partial single-joint 
bends. 
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Figure 6: A polygonal chain that cannot be opened with single-joint moves. 



Proof: The basic idea is similar to the one in Theorem 3.1 and also estab- 
lishes a reduction of Partition. (Refer to Figure for an overview.) As before, 
we write bi for the joints, and Si = bi) for the segments. We use the idea 

of the construction from Lemma |3.3| to construct a rigid frame, with the key 
corresponding to the free end of that chain. The frame has one end, bo, of the 
polygonal chain wedged into the corner &13, which has angle ip <^ tt/2. Because 
of the degeneracy at &13, none of the joints &i, . . . , 612 can be moved individually 
without causing a self-intersection between bo and the chain in the neighbor- 
hood of 613: Just like in the proof of Lemma 3.3, none of the joints &i, . . . , 612 
lies in the area of possible locations of feasible rotations. This continues to be 
the case while (j) -\- < t: , i.e., while the sum of angles at 613 and at 617 does 
not change significantly. 

Again, the "key" contains the n segments Sig, . . . , Sig+n of integral lengths 
oi, . . . , a„ that encode an instance of Partition. As before, let S denote the set 
of integers for the Partition instance. We also use "long" auxiliary segments 
of lengths L/2 and L, where L 3> = Here segments sig and S2o+n 

have length L/2, while S2i+n has length L. 

The critical dimensions of the frame are chosen such that the key can just 
be removed from the frame if and only if it can be collapsed to a length of L. 
Removing the key consists in pulling it through the narrow bottleneck formed 
by the segments S3 and sg by extending the "spring" formed by S14 and S15, 
while moving the "keyholder" sis down by a distance of L -I- e. This is possible 
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if and only if there is a feasible partition. 



(degeneracy) 




angle 9 



[ key j 
(vertical size L ) 



( spring) 

(segments of length 
(U2cos^) + 0(e)) 




angle V b^^ 



(bottleneck) 



Figure 7: Illustration of the proof of Theorem 3.4. Note that lengths are not 
drawn to scale, in order to show sufficient details; in particular, the dimensions of 
the bottleneck are much smaller than the edges encoding the partition instance. 



More precisely, let the angle at 617 be -0 <C tt. Let e — 0{l/nL). We assume 
that the dimensions of the frame are chosen sufficiently large to guarantee that 
moving bn down by a vertical distance of L + £ increases its distance from 
613 by Lcosip + &{e), i.e., the angles at 617 and at 613 do not change much. 
The segments su and S15 forming the spring have length L/2cosV' + 6(e), so 
extending the spring will just suffice to move the keyholder sis down by L + e, 
but not more. The vertical "height" of the bottleneck, i.e., the length of the 
segments S3 and sg, is e/3, while the horizontal "width" xg — 2:3 = xg — X2 — 
is significantly smaller. (As we will discuss below, this forces the keyholder to be 
roughly vertical throughout the motion.) For the initial position of the key inside 
of the frame, we assume L < j/ig — j/3 < L + e/3 and L < y2o+n —ys < L + e/3. 
Finally, ye ~ yr — L + ld{e) and xq — X5 = 6(^), so the dimensions of the 
rectangle formed by &4,65,&6j&7 are not large enough to change the basically 
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vertical orientation of the location segments sig, si9+„, and S2o+n- 

Now assume that there is a feasible partition, S = Si U S2, such that 
X^ieSi '^i — 12i£S2 performing a (finite) "wiggly" sequence of moves, we 

can move the partition segments such that (1) any segment sig+i representing 
tti G Si satisfies yi — yi-i — ai + 0{e^) and \xi — Xi-i\ — O(e^), so that sig+i 
is pointing up; (2) any segment sig+i representing € S2 satisfies yi^i ~ yi = 
tti + 0{e^) and \xi — Xi-i\ = O(e^), so that sig+i is pointing down; and (3) 
we end up placing 520+n within Euclidean distance 0{e^) from big and placing 
621+n at distance e/3 + 0{e^) from 63 and bg,. Thus, extending the spring by 
an appropriate wiggly motion moves the key through the bottleneck. Now it is 
easy to open the joint 613, and unfold the whole chain. 

Conversely, assume that the chain can be unfolded. As discussed above, the 
sum of angles at 613 and at bn has to change significantly before the frame 
ceases to be rigid. Now note that the dimensions of the bottleneck force the 
keyholder segment to be roughly vertical, i.e., to have slope within 0{l/e^) of 
vertical. (See Figure^.) Furthermore, we noted above that any feasible vertical 
motion of 617 does not change the angle at 617 by a significant amount; it is 
clear that this also prevents the angle at 613 from changing much. Therefore, 
the frame remains rigid until the key has been removed from the lock. 



0(Le3) 



(entering 
bottleneck) 



keyholder s^^ 




(leaving 
bottleneck) 



0(Le2) 



Figure 8: All segments have to be inside of the shaded region when moving 
through the bottleneck, i.e., must be close to being vertical. (Horizontal scale 
and size of the bottleneck are vastly exaggerated to allow sufficient resolution. 
In scale, the shaded region is basically a vertical line.) 
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Now consider the positions of segments sis and S2i+n when 52i+n crosses 
the horizontal hne y = y2- By the dimensions of the bottleneck, S2i+n must 
have a slope within r2(l/e'^) of vertical. Furthermore, by construction of the 
rectangle 646566^7, we are assured that yis — y20+n < 0(e), i.e., S20+n cannot be 
significantly below sis- On the other hand, 620+4^ must be below the horizontal 
line y = ya when bis has been moved down by a vertical distance of L + e. 
Since no segment within the key can change its vertical slope significantly while 
S21+TI is within the bottleneck (they must all remain wedged between sig and 
S21+TI and be strictly contained in the shaded region in Figure we conclude 
that 2/20+n — 2/18 < 0{e) upon leaving the bottleneck, i.e., S20+n cannot be 
significantly above sis- 

Therefore, the sets 6*1 = {i G 1, . . . , n\yi > yi-i} and ^2 = {« G 1, . . . , n\yi < 
yi-i} upon entering the bottleneck from above and leaving it from below must 
satisfy J2ieSi = ^ieS2 Since e ^ 1, this implies that there is a 

feasible partition. □ 

For the case of monotonic bend operations, the above proof can be easily 
modified: 

Corollary 3.5 It is NP-hard to decide if a polygonal chain P with a single 
vertex-to-vertex incidence can he straightened by monotonic partial single-joint 
bends. 

Proof: The joints in the construction shown in Figure ^ that may not be 
changed monotonically are &14, . . . , 620+n: the ones that are not part of the 
frame. By using small quadruple joint gadgets as in the construction for Theo- 
rem we get a chain that can be opened with monotonic moves, if and only 
if it can be opened. □ 
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4 Algorithms 



In this section, we turn our attention to positive algorithmic results, giving 
efBcient algorithms for deciding if particular bend sequences are feasible. We 
consider here only the case of complete bends. 

Consider an arbitrary permutation, a — of the bends along 

a wire. In order for cr to be a foldable sequence, it is necessary and suffi- 
cient that for each j ~ 1, . . . ,n the bend hi- is foldable. Recall that in our 
notation P({ii, . . . denotes the partially bent chain after the bends at 

^ii, ■ ■ • have been straightened. The point hi. splits P{{i\, . . . into 

two subchains; let Pq (resp., Pn+i) denote the subchain containing the end- 
point bo (resp., &„+i). Now, hi^ is foldable if the joint at hi. can be straightened 
without causing a collision to occur between Pq and Pn+i at any time during 
the rotation about hi. . We can assume, without loss of generality, that Pq is 
fixed and that Pn+i is pivoted about hi^ . During this bend operation at hi^ , 
each point, u, on Pn+i moves along a circular arc, A^, subtending an angle 9i ., 
centered on hi-. It is clear that in order for the bend to be feasible, none of 
these arcs Au may cross the chain Pq, for all choices of points u on Pn+i- 

If the perpendicular projection of hi . onto the line containing an edge e of 
Pn+i lies on the edge, let 6 e denote the projection point. (Each edge of 
Pn+i has at most one projection point.) Let U denote the union of the set of 
vertices of Pn+i and the set of projection points on edges of Pn+i- In the lemma 
below, we observe that, in order to test feasibility of straightening the bend hi . , 
it suffices to consider only the feasibility of the final position of the chain Pn+i 
and to test Pq for intersection with the discrete set A ~ {Au : u e U}. See 
Figure ^. 

Lemma 4.1 Joint hi. is foldable if and only if (1) no arc of A intersects Pq, 
and (2) after the bend, no segment of Pn+i intersects a segment of Pq- 

Proof: If joint bi is foldable then, by definition, there can be no intersection 
of Pn+i with Pq during its rotation about hi,. This implies conditions (1) and 
(2). 

If conditions (1) and (2) hold, then we claim that there can be no intersection 
of Pn+i with Pq during the rotation. Consider a subsegment, s, of Pn+i whose 
endpoints are consecutive points of U. (Note that at least one endpoint of s 
must be a vertex of P„+i.) During the rotation, it sweeps a region Rg that is 
bounded by two circular arcs centered at bi^ , corresponding to the trajectories 
of its endpoints during the rotation, and two line segments, corresponding to 
the positions of s before and after the rotation. Here, we are using the fact that 
the distance from bi^ to a point q E s monotonically changes as a function of 
the position of q on s. (The projection points were introduced in order to assure 
this property.) Our claim follows from the fact that Pq is a simple, connected 
chain: It cannot intersect s at some intermediate stage of the rotation unless it 
intersects the boundary of the region Rg. Such an intersection is exactly what 
is being checked with conditions (1) and (2). □ 
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Figure 9: Foldability of the joint 6^: The subchain P„+i is shown with thicker 
lines (two dashed copies show it after different stages of rotation about bi.). 
Each vertex and each projection point (shown as black disks) of Pn+i moves 
along a circular arc, shown using a thin dashed arc. In this example, the rotation 
shown is not feasible, as it fails both conditions (1) and (2) of the lemma. 



Lemma 4.2 For any S Q B, and any hi ^ S, one can decide in O(nlogn) time 
if joint bi is foldable for the chain P{S). 

Proof: Using standard plane sweep methods for segment intersections, 
adapted to include circular arcs, we can check in O(nlogn) time both con- 



ditions ((1) and (2)) of Lemma 4.1. Events in the sweep algorithm correspond 
to joints and to vertical points of tangency of circular arcs, assuming we use 
a vertical sweep line. During the sweep, we keep track of the vertical ordering 
of the segments and arcs that cross the sweep line; we check for intersection 
between any two objects that become adjacent in this ordering, stopping if a 
crossing is detected. Since we process 0{n) events, each at a cost of 0{\ogn), 
the time bound follows. □ 

Remark 4.3 Condition (2) can be tested in linear time, by Chazelle's trian- 
gulation algorithm. We suspect that condition (1) can also be tested in linear 
time. Condition (1) involves testing for rotational separability of two simple 
chains about a fixed center point (bi ), which is essentially a polar coordinate 
variant of translational separability ( which is easily tested for simple chains us- 
ing linear-time visibility (lower envelope) calculation). The issue that must be 
addressed for our problem, though, is the "wrap-around" effect of the rotation; 
we believe that this can be resolved and that this idea should lead to a reduction 
in running time of a factor oflogn. 

Corollary 4.4 The foldability of a permutation a can be tested in 0{n^ ^ogn) 
time. 
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We obtain improved time bounds for testing the feasibility of a particularly 
important folding sequence: the identity permutation. Many real tube-bending 
and wire-bending machines operate in this way, making bends sequentially along 
the wire/tube. (Such is the case for the hydraulic tube-bending machines at 
Boeing's factory, where this problem was first suggested to us.) Of course, there 
are chains P that can be straightened using an appropriate folding sequence 
but cannot be straightened using an identity permutation folding sequence; see 
Figure |^(e). However, for this special case of identity permutations, we obtain 
an algorithm for determining feasibility that runs in nearly linear time: 

Theorem 4.5 In time O(nlog^n) one can verify if the identity permutation 
(a = (1,2,.. .,n) J is a foldahle permutation for P. 

Proof: For notational convenience, we consider the equivalent problem of 
verifying if it is feasible, in the order &i, 62, . . . , to bend the joints bi from 
joint angle tt to final angle 0i, thereby transforming a straight wire into the final 
shape P, rather than our convention until now of considering the problem of 
performing bend operations to straighten the chain P. 

Thus, consider performing the bends in the order given by the identity per- 
mutation a, and consider the moment when we are testing the foldability of bi. 
The subchain Pn+i, from bi to is a single line segment, bib„+i, since the 

joints . . . ,bn are straight at the moment. Thus, verifying the foldability 
of bend bi amounts to testing if the segment bibn+i can be rotated about bi 
by the desired amount, without colliding with any other parts of the subchain 
Pq of P, from 60 to bi. In other words, we must do a wedge emptiness query 
with respect to Pq, defined by bi, segment bibn+i, and the angle 9i. Since Pq is 
connected, emptiness can be tested by verifying that the boundary of the wedge 
does not intersect Pq. (See Figure |l^.) Thus, we can perform this query by 
using (straight) ray shooting and circular-arc ray shooting in Pq; the important 
issue is that Pq is dynamically changing as we proceed with more bends. 

However, in order to avoid the development of potentially complex dynamic 
circular-arc ray shooting data structures, we devise a simple and efficient method 
that "walks" along portions of Pq, testing for intersection with the circular arc, 
7, from bn+i to 6Ji_|_i, where 6^+1 is the location of 6,1+1 after the bend at bi has 
been performed. 

In particular, we keep track of a "painted" portion of Pq, which corresponds 
to the subset of Pq that has been "walked over." We consider the chain Pq to be 
a degenerate simple polygon, having two sides which form a counterclockwise 
loop around Pq. We consider the case in which the bend at bi is a rotation of 
the segment clockwise to the segment bib'^^^i, the case of a counterclock- 

wise bend at bi is handled similarly. When we perform a bend at bi, we walk 
(counterclockwise) along the unpainted portions of Pq, between two points, a 
and a', on the boundary of Pq, where a and a' are defined according to cases 
that depend on the outcomes of two ray-shooting queries: 

(a) If both of the rays bibn+i and &i6'„+i miss Pq (and go off to infinity), then 
there is nothing more to check: the rotation at bi can be done without 
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Figure 10: Testing the foldability of the joint 6^. This example is intended to 
illustrate a generic step in the algorithm; for this particular chain, note that it 
is not feasible to make the bends bi, . . . ,bi-i to get to the state shown. 




Figure 11: Case (a): Both of the rays bibn+i and bib'^_^_i miss Pq and go off to 
infinity. 
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interference with Pq, since Po is a (connected) polygonal chain lying in the 
complement of the wedge defined by bibn+i and bih'^j^i. See Figure [Tl| . 




Figure 12: Case (b): Both of the rays bibn+i and bib'^+i hit Po- The walk 
extends from a to a' over the highlighted portion of Pq, painting any previously 
unpainted portion of it. 



(b) If both of the rays bibn+i and hit Pq, then we let a and a' (respec- 

tively) be the points on the boundary of Pq where they first hit Po. See 
Figure |l^. 

(c) If exactly one of the rays bibn+i and foj&Ii+i hits Pq while the other misses 

Po (and goes off to infinity), then we define a and a' as follows. Assume 

that the ray bib^+i hits Po (and the ray misses Pq); the other case 

is handled similarly (see Figure |l^, right). Then, we define a to be the 
point on the boundary of Pq where the ray bibn+i hits Pq, and we define 
a' to be the point on the boundary of Po where a ray from infinity in the 
direction 6„+ifei (towards bi) hits Pq. See Figure [T^, left. 

During the walk from a to a' along the boundary of Po, we test each segment 
for intersection with the circular arc 7 in time 0(1). Whenever we reach a 
portion of the boundary that is already painted, we skip over that portion, 
going immediately to its end. Already painted portions have endpoints that 
were determined by rays in previous steps of the painting procedure. Since 
there are only a total of 0{n) rays (one per edge of P), this implies only 0{n) 
endpoints of painted portions. As we walk, we mark the corresponding portions 
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Figure 13: Case (c): Exactly one of the rays bibn+i and hits Pq: bibn+i 

hits Pq (left) or bib'^_^_i hits Pq (right). The walk extends from a to a' over the 
highlighted portion of Pq, painting any previously unpainted portion of it. 



over which we walk as "painted." Since, by continuity, it is easy to see that 
the painted portion of any one segment of Pq is connected, we know that we 
must encounter at least one vertex of Po between the time that the walk leaves 
a painted portion and the time that the walk enters the next painted portion. 
Thus, during a walk, we charge the tests that we do for intersection with 7 off 
to the vertices that are being painted. The remainder of the justification of the 
algorithm is based on two simple claims: 

Claim 4.6 There is no need to walk back over a painted portion in order to 
check for intersections with an arc 7 at some later stage. 

Proof: The fact that we need not walk over a painted portion testing again 
for intersections with 7 follows from the fact that with each bend in the sequence, 
the length of the segment bibn+i that wc arc rotating goes down by the length 
of the last link. Thus, if the motion of the tip, bn+i, sweeps an arc 7 that does 
not reach a portion /j, of the boundary of Pq when the link bibn+i is straight, 
it cannot later be that a link bjbn+i {j > i) can permit the tip bn+i to reach 
the same portion fj, when pivoting is done about bj; this is a consequence of the 
triangle inequality. □ 

Claim 4.7 In testing for intersection with 7, we check enough of the chain Pq." 
if any part of it intersects 7, then it must lie on the portion between a and a' 
over which we walk. 

Proof: In case (a), there is nothing to check. In case (b), the closed Jordan 
curve from bi to a (along a straight segment), then along the boundary of 
the simple polygon Pq to a', then back to bi (along a straight segment) forms 
the boundary of a region whose only intersection with Pg is along the shared 
boundary from a to a'; thus, if 7 lies within this region (i.e., does not intersect 
the boundary of Pq from a to a'), then 7 does not intersect any other portion 
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of Po- In case (c) we argue similarly, but we use the Jordan region defined by 
the segment from bi to a, the boundary of Pq from a to a', the ray from a' to 
infinity (in the direction of bibn+i), then the reverse of the ray 6i6^_|_]^ back to bi. 
□ 

The total time for walking along the chain Pq can be charged off to the 
vertices of P, resulting in time 0{n) for tests of intersection with arcs 7, exclusive 
of the ray shooting time. The final time bound is then dominated by the time 
to perform n straight ray shooting queries in a dynamic data structure for the 
changing polygonal chain Pq] these ray shooting queries are utilized both in 
testing for intersection with the segment ^i^Jj+i and in determining the points 
a and a' that define the walk. These ray-shooting queries and updates are done 
in time 0(log^ n) each, using existing techniques (||l^), leading to the claimed 
overall time bound. □ 

Next, we turn to two other important classes of permutations. Again, for 
notational convenience, we consider the problem of verifying if it is feasible, 
in the order given by the permutation, to bend the joints bi from joint angle 
TT to final angle 9i, thereby transforming a straight wire into the final shape 
P. We say that a permutation is an outwards folding sequence (resp., inwards 
folding sequence) if at any stage of the folding, the set of bends that have been 
completed, and therefore are not straight, is a subinterval, bi, bi+i, . . . ,bj (resp., 
a pair of intervals &i, 62, • ■ ■ , and bj, . . . , bn); thus, the next bend to be 
performed is either bi-i or (resp., fe^+i or Inwards and outwards 

folding sequences are a subclass of permutations that model a constraint im- 
posed by some forming machines. See Figure The identity permutation is 
a folding sequence that is a special case of both an inwards and an outwards 
folding sequence. 

We show that one can efficiently search for a folding sequence that is inwards 
or outwards. Our algorithms are based on dynamic programming. 

First, consider the case of outwards folding sequences. We keep track of 
the state as the pair (z,j) representing the interval of bends (&i, fo^+i, . . . , 6j) 
already completed. We construct a graph Q whose O(n^) nodes are the states 
(i, j) (with I < i < j < n) and whose edges link states that correspond to 
the action of completing a bend at bi-i or bj+i (if the starting state is {i,j))- 
Thus, each node has constant degree. Our goal is to determine if there is a 
path in this graph from some (i, i), for i G {1, 2, . . . , n} to (1, n). Wc augment 
this graph with a special node vq, linked to each node (i, i). Then, our problem 
is readily solved in O(n^) time once we have the graph constructed, since it is 
simply searching for a path from node i^o to node (l,n). (Alternatively, we can 
construct the graph as we search the graph for a path.) In order to construct 
the graph, we need to test whether bend or bj^i can be performed without 
intersecting the folded chain, P' , linking to 6j+i. This is done in a manner 
very similar to that we described above for the case of identity permutations: 
we perform ray-shooting queries in time 0(log^ n) and then use a "painting" 
procedure to keep track of the states of 2n — 2 "walks" that determine circular- 
arc ray shooting queries. In particular, there is a separate painting procedure 
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Figure 14: An intermediate state in the bending of an outwards (left) and 
an inwards (right) folding sequence. For the outwards folding sequence on the 
left, the next bend is either 6i„i or the new positions of the chain are 

shown dashed. For the inwards folding sequence on the right, the next bend is 
either bi+i or 



corresponding to each of the n — 1 choices of i and to each of the n — 1 choices of 
j. For example, for a fixed choice of i, the painting procedure will consider each 
of the possible bends . . . ,bn in order, allowing us to amortize the cost of 
checking for intersections with the circular arc 7 associated with each bend. In 
total, the cost of the walks is 0{n^), while there may also be O(n^) ray shooting 
queries (in a dynamically changing polygon). Thus, the total cost is dominated 
by the ray shooting queries, giving an overall time bound of 0(n^ logn). 

For the case of an inwards folding sequence, we build a similar state graph 
and search it. However, the cost of testing if a bend is feasible is somewhat 
higher, as we do not have an especially efficient procedure for testing the fold- 
ability of a polygonal chain. (Our painting procedures exploit the fact that the 
link being folded is straight.) Thus, we apply the relatively naive method of 



testing feasibility given in Lemma 4.2, at a cost of 0{nlogn) per test (which 



potentially improves to 0{n) time, if our conjecture mentioned in the remark 
after the Lemma is true). Thus, the overall cost of the algorithm is dominated 
by the O(n^) feasibility tests, at a total cost of 0{n^ \ogn). In smumary, we 
have: 

Theorem 4.8 In time 0{n^ log^ n) one can determine if there is an outwards 
folding sequence; in time 0(n^ logn) one can determine if there is an inwards 
folding sequence. 
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5 Conclusion 



We conclude with some open problems that are suggested by our work: 

(1) Is the bend sequencing problem for wire folding strongly NP-complete, or 

is there a pseudo-polynomial-time algorithm? If not in wire bending, is it 
strongly NP-complete for the 3-dimensional sheet metal folding problem? 

(2) Is it NP-hard to decide if a polygonal chain in three dimensions can be 

straightened? In |^ simple examples of locked chains in three dimensions 
are shown; can these be extended to a hardness proof for the decision 
problem? 

(3) In practice, in order to make a bend using a punch and die on a press 

brake, it is necessary to consider accessibility constraints. For each bend 
operation, the die is placed on one side of the material, while the punch is 
placed on the other side. The bend is formed by pushing the punch into 
the die (which has a matching shape) , with the material in between. (See 
Wang [|^.) In the simplest model of this operation on a wire, we can 
consider the punch and the die to be oppositely directed rays that form 
a bend by coming together (from opposite sides of the wire) so that their 
apices meet at the bend point. The accessibility constraint in this simple 
model is that the rays representing the punch and die must be disjoint 
from the wire structure both at the initial placement of these "tools" and 
during the bend operation itself. 

(4) Can the foldability of a permutation be decided in subquadratic time for 

wire bending? This would be possible if one had a dynamic data structure 
that will permit efhcient (sublinear) queries for the foldability of a vertex. 
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